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With a sample of about 384 million BB pairs recorded with the BABAR detector, we search for the flavor- 
changing charged current transition B^ — > r^f and perform an amplitude analysis of the effective flavor- 
changing neutral current transition B^ ip{1020)K* (892)^ . We also extend our search for other K* final 
states in the decay — ► ip{1020)K''^ with a large K*'^ — > K^n~ invariant mass. Two samples of events with 
one reconstructed hadronic B decay or one reconstructed semileptonic B decay are selected, and in the recoil 
a search for B^ — > t^u is performed. We find a 2.6 tr (3.2 a not including expected background uncertainty) 
excess in data which can be converted to a preliminary branching fraction central value of B{B^ t^u) = 
(1.20lo;38lo 3Q ± 0.22) X 10-". With the decay _B± v'(1020)/s:* (892)± , twelve parameters are measured, 
where our measurements of /j, = 0.49 ± 0.05 ± 0.03, f± = 0.21 ± 0.05 ± 0.02, and the strong phases point to 
the presence of a substantial helicity-plus amplitude from a presently unknown source. 
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Until the new frontier energy is available at the LHC 
collider, the Standard Model is being investigated in 
detail via the rich flavor structure of the quarks and 
leptons. The study of rare processes, such as flavor- 
changing neutral and charged currents, already pro- 
vides constraints on new physics and will allow us 
to disentangle the flavor structure of particles at the 
new energy scale once LHC energy is accessible. Until 
2008, flavor physics of the third-generation particles, 
such as &-quark, offers best prospects for discoveries. 
For example, virtual "loop" transitions of the B me- 
son involve the heaviest presently known particles in 
the FCNC loop and may reveal new particles through 
subtle effects. In this paper, we report the most recent 
("hot") results on two rare processes B^ t^v and 
B^ -> ip{1020)K*{892)^ with the BaBAR detector. 

In the Standard Model (SM), the purely leptonic de- 
cay B^ T^v proceeds via quark annihilation into a 
boson, as shown in Fig. 1. The branching fraction 



is given by: 



B{B^ 



Stt 



(1) 

where Vu,b is a quark mixing matrix element [1, 2], 
/b is the B meson decay constant, Gi;- is the Fermi 
constant, is the B~^ lifetime, and and m-r are 
the 5+ meson and r masses. Physics beyond the SM, 
such as two-Higgs doublet models, could enhance or 
suppress B{B'^ — > r^Vr) through the introduction of 
a charged Higgs boson [3-5], see Fig. 2. Branching 
fraction in Eq. (1) is modified by a factor: 



1 - tan^ /3 



''H+ 



(2) 



where tan (3 is the ratio of vacuum expectation values 
of the two Higgs doublets. Using theoretical calcula- 
tions of /b from lattice QCD and experimental mea- 
surements of \Vub\ from semileptonic B decays, this 
purely leptonic S-decay can be used to constrain the 





Figure 1; Feynman diagram describing the 
decay. 



Figure 2: Feynman diagram with potential charged Higgs 
contribution to the — > t^Vt decay. 
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parameters of theories beyond the SM. Or, assuming 
that SM processes dominate and using the value of 
\Vub\ determined from semileptonic _B-decays, purely 
leptonic decays provide a clean experimental method 
of measuring fs precisely. 

The branching fractions for _B+ fi^iy and 
e'^v are suppressed by factors of ^ 5 x 10^^ and 
- 10""^ with respect to B+ ^ r+Vr- The SM es- 
timate of the branching fraction for B^ t^i^t, 
using iKbl = (4.31 ± 0.30) x lO-^ [6] and fs = 
0.216 ± 0.022 GeV [7] in Eq. 1 is (1.6 ± 0.4) x lO""*. 
However, a search for _B+ t^Vt is experimen- 
tally challenging due to the large missing momen- 
tum from multiple neutrinos, which makes the signa- 
ture less distinctive than in the other leptonic modes. 
In a previously published analysis using a sample 
of 223 X 10^ T{AS) decays, the BaBAR collabora- 
tion set an upper limit of B{B^ T^Vr) < 1.8 x 
10"'' at the 90% confidence level (CL) [8]. The Belle 
Collaboration has reported evidence from a search for 
this decay where the branching fraction was measured 
tohe B{B+ T+Vr) = (1.79t°i^(stat.)t°-^^(syst))x 
10-4 [9]. 

A vector- vector _B-meson decay, such as -B ^ ipK* , 
is characterized by three complex helicity amplitudes 
Ai), which correspond to helicity states A = —1, 0, +1 
of the vector mesons. The Aiq amplitude is ex- 
pected to dominate [10] due to the {V — A) nature 
of the weak interactions and helicity conservation in 
the strong interactions, see Fig. 3. A large fraction of 
transverse polarization observed by BaBAR and con- 
firmed by Belle, along with more recent measurements 
of polarization in rare vector-vector B meson decays 
B (j)K* and pK*, indicate a significant departure 
from the expected predominance of the longitudinal 
amplitude [11-15]. The rate, polarization, and CP 
measurements of B meson decays to particles with 
nonzero spin are sensitive to both strong and weak 
interaction dynamics and are discussed in a recent re- 
view [6, 16]. 

The polarization anomaly in vector-vector B me- 
son decays suggests other contributions to the decay 
amplitude, previously neglected. This has motivated 
a number of proposed contributions from physics be- 
yond the standard model [17]. Depending on New 
Physics model, hierarchy of decay amplitudes could 
be modified, as shown in Fig. 4. In addition, there 
are new mechanisms within the standard model which 
have been proposed to address the anomaly, such as 
annihilation penguin [18] or electroweak penguin, or 
QCD rescattering [19], 

In order to distinguish the models, the BaBAR ex- 
periment extended the study of the i?" (j)K*^ de- 
cays with the tensor (J^ = 2+), vector (J^ = 1~), 
and scalar (J^ = 0+) K*° [15]. The vector-tensor 
results are in agreement with quark spin-flip suppres- 
sion [10] and ^0 amplitude dominance, whereas the 
vector- vector mode contains substantial ampli- 




Figure 3: Feynman diagram describing the B tpK* de- 
cay. Due to the {V — A) nature of the weak interactions, 
helicity conservation in the strong interactions, and quark- 
spin suppression shown in the diagram, we expect the hi- 
erarchy Ao > A+i > A-i. 




Figure 4: Feynman diagram describing the B ipK* 
decay. Scalar interaction may modify the hierarchy to 
A+i S> ^0 ^ ^-1, while supersymmetric interaction with 
(V + A) couplings could produce Aq S> A^i ^ ^l+i. 



tude, corresponding to anomalously large transverse 
polarization. 

We now investigate the polarization puzzle with a 
fuU amphtude analysis of the ipK* (892)^ de- 

cay. In this paper, we report twelve independent pa- 
rameters for the three B'^ and three B~ decay am- 
plitudes, six of which are presented for the first time. 
Moreover, we use the dependence on the Kt: invari- 
ant mass of the interference between the = 1^ and 
0+ {Ktt)^ components [15, 20, 21] to resolve the dis- 
crete ambiguity between the Ai^i and Ai^i helicity 
amplitudes. 
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Figure 5: Diagram of he —* r^u reconstruction tech- 
nique. 



2. The BaBar Detector 

We use a sample of 383.6±4.2 million T(4S') BB 
events collected with the BaBAR detector [22] at the 
PEP-II e~^e~ asymmetric-energy storage rings. The 
e+e~ center-of-mass energy y/s is equal to 10.58 GeV. 
The sample corresponds to an integrated luminosity 
of 346fb~^ at the T(45') resonance (on-resonance) and 
36.3 fb~^ taken at 40MeV below the BB production 
threshold (off-resonance) which is used to study back- 
ground from e+e~ ^ // (/ = u, d, s, c, r) continuum 
events. The detector components used in this analysis 
are the tracking system composed of a five-layer silicon 
vertex detector and a 40-layer drift chamber (DCH), 
the Cherenkov detector for charged tt-K discrimina- 
tion, a Csl calorimeter (EMC) for photon and elec- 
tron identification, and an 18-layer flux return (IFR) 
located outside of the 1.5 T solenoidal coil and instru- 
mented with resistive plate chambers for muon and 
neutral hadron identification. For the most recent 133 
fb~^ of data, a portion of the resistive plate chambers 
has been replaced with limited streamer tubes. 

A GEANT4-based [23] Monte Carlo (MC) simula- 
tion is used to model signal efficiencies and physics 
backgrounds. The r lepton decay is modeled using 
EvtGen [24]. Beam-related background and detector 
noise from data are overlaid on the simulated events. 
Simulation samples equivalent to approximately three 
times the accumulated data are used to model BB 
events, and samples equivalent to approximately 1.5 
times the accumulated data are used to model con- 
tinuum events. We determine selection efficiencies for 
signal events using a MC simulation where one B~^ 
meson decays to r+z/, while the other is allowed to 
decay into any final state. 

3. B± T^v 

First we discuss the new preliminary analysis of the 
decay. Due to the presence of multiple 



neutrinos, the i?+ —t t^Vt decay mode lacks the kine- 
matic constraints that are usually exploited in B de- 
cay searches in order to reject both continuum and 
BB backgrounds. The strategy adopted for this anal- 
ysis is to reconstruct exclusively the decay of one of 
the B mesons in the event, referred to as the "tag" 
B. The remaining particle(s) in the event (the "re- 
coiling system" ) are assumed to come from the other 
B and are compared with the signature expected for 
B^ t'^Vt, see Fig. 5. There are two "tag" methods 
used in this analysis: full reconstruction of a hadronic 
B meson decay final state or partial reconstruction of 
a semileptonic decay. We will discuss in detail the new 
results with the former method and then briefly show 
updated results of the latter method with slightly im- 
proved precision compared to last year. Then we will 
combine the two results. 

3.1. Analysis Method 

In order to avoid experimenter bias, the signal re- 
gion in data is blinded until the final yield extraction 
is performed. The B'^ — > t'^Vt signal is searched for 
in both leptonic and hadronic t decay modes consti- 
tuting approximately 71% of the total t decay width: 
r"*" e'^vv, T"*" — > fi^iyV, t+ —>■ 7r+77, and t+ —^ 
tt^tt'^V. We do not consider the r+ tt+tt^tt+F 
mode since we found it to be dominated by back- 
ground events. 

The tag B candidate is reconstructed in the set 
of hadronic B decay modes B~ D'^*^^X^ [25], 
where X~ denotes a system of charged and neu- 
tral hadrons with total charge —1, composed of 
UiTT^ n2K^ n^K^ n^TT^ , where ni + n2 < 5, 71.3 < 2, 
and 714 ^ 2. 

The selected sample of tag B candidates is used 
as normalization for the determination of the branch- 
ing fraction. Wc reconstruct D*'^ D^tt^ , D^j; 
D° -> K-TT+,K-n+TT^, X-7r+7r-7r+, K^TT+n-; and 
TT+TT^. The kinematic consistency of tag B 
candidates is checked with two variables, the beam 
energy-substituted mass toes = ■\/s/4 — p| and the 
energy difference AE = Eb — ^fs/"^- Here ^/s is 
the total energy in the T{AS) center-of-mass (CM) 
frame, and pb and Eb denote, respectively, the mo- 
mentum and energy of the tag B candidate in the 
CM frame. The resolution on AE is measured to 
be (TAB = 10 — 35MeV, depending on the decay 
mode, and we require \AE\ < 3aAE- For each recon- 
structed B decay mode, its purity V was estimated as 
the ratio of the number of peaking events with mss 
> 5.27GeV/c to the total number of events in the 
same range, and is evaluated on on-resonance data. 
In the event of multiple tag B candidates being recon- 
structed, the one with the best purity V is selected. 

The background consists of e'^e^ ^ il ~ 
u,d,s,c) events and of other T{4S) B^B^ or 



fpcp07_lll 



4 



Flavor Physics and CP Violation Conference, Bled, 2007 




beam axis and less than 1.5 cm transverse to it. 
The T lepton is identified in the four decay modes: 



5.28 5.29 5.3 
m^s (GeV/c^) 

Figure 6: Cumulative distribution of the energy substi- 
tuted mass rriES of the tag B candidates in data. The non 
peaking background components are added one on top of 
the other. From top to bottom are B~ , B^B'^ , uds and 
cc. 



B^B^ decays, in which the tag B candidate is mistak- 
enly reconstructed from particles coming fr'om both 
B mesons in the event. To significantly reduce the 
e"'"e~ — > qq background we require the angle de- 
fined in the CM frame, between the thrust axis [26] 
of the tag B candidate and the thrust axis of all the 
charged and neutral reconstructed candidates in the 
event excluding the ones that form the to satisfy 
the requirement jcos^^^^l < 0.9. 

In order to determine the number of correctly re- 
constructed B'^ we classify the background events in 
four categories: e+e~ cc; e+e~ uu, dd, ss; 
B'^B^; and B~^B^ . The shapes of these background 
distributions are taken from MC simulation. The nor- 
malizations of the e"'"e~ —^ cc and e~^e~ — *■ uu, dd, ss 
backgrounds are taken from off-resonance data, scaled 
by the luminosity. The normalization of the B^B^ , 
B^B^ components are instead obtained by means of 
a fit to the tties distribution in a sideband re- 
gion (5.22GeV/c2 < mEs < 5.26GeV/c2). The back- 
ground contamination in the signal region (tties > 
5.27GeV/c^) is extrapolated from the fit and sub- 
tracted from the data. We estimate the total number 
of tagged B's in the data as Nb = (5.92 ± O.llstat.) x 
10^. Fig. 6 shows the tag B candidate tties distri- 
bution in data compared with MC, and the obtained 
background subtracted mEs distributions in data. 

After the reconstruction of the tag B meson, a se- 
lection is applied on the recoiling system in order to 
enhance the sensitivity to B^ — > t^Vt decays. We 
require the presence of a well reconstructed charged 
track (signal track) with electric charge opposite to 
the tag B. The signal track is required to have at 
least 12 hits in the DCH, its momentum transverse 
to the beam axis, p^, is required to be greater than 
0.1 GeV/c, and its point of closest approach to the in- 
teraction point must be less than 10.0 cm along the 



r+ 



r + 



T' e' vv, T' fi^vV, r+ — )■ 7r+I7, and 
7r+7r'^77. Particle identification criteria on the signal 
track are used to separate the four categories. The 
T+ — > 7r+7r"77 sample is obtained by associating the 
signal track, identified as pion, to a tt" reconstructed 
from a pair of neutral clusters with invariant mass 
between 0.115 and 0.155 GeV/c^ and total energy 
greater than 250 MeV. In case of multiple 7r+7r° candi- 
dates the one with largest momentum p*+^o is chosen. 

Most of the background due to continuum events 
and incorrectly reconstructed tag B candidates (com- 
binatorial) is rejected by requiring a mode dependent 
cut on \cos9^g\. Most of the remaining sources of 
background consists of B'^B~ events in which the 
tag B meson was correctly reconstructed and the re- 
coil contains one track and additional particles which 
are not reconstructed by the tracking detectors and 
calorimeter. From MC simulation we observe that 
most of this background is from semileptonic B de- 
cays. 

We define the discriminating variable -Boxtra as the 
sum of the energies of the neutral clusters not asso- 
ciated with the tag B or with the signal tt" of the 
r+ — > tt+tt'^F mode. For neutral clusters contributing 
to E'cxtra we need to determine mode-dependent min- 
imum energy thresholds, shown in Table I, in order to 
ensure a good data-MC agreement at the preselection 
level, which is a looser selection using only particle 
identification criteria and a moderate cut on the can- 
didate's momentum. Signal events tend to peak at 
low i?oxtra values whereas background events, which 
contain additional sources of neutral clusters, are dis- 
tributed toward higher £^cxtra values. 



Table I EMC barrel and EMC endcap minimum cluster 
energies ( MeV) for electrons, muons and pions samples 



Sample 


barrel 


endcap 


electrons 


65 


70 


muons 


50 


55 


pions 


50 


70 



Other variables used to discriminate between signal 
and background are the CM momentum of the signal 
candidates, the multiplicity of extra charged track(s) 
and 7r° (s) in the recoil and the direction of the miss- 
ing momentum four-vector in the CM frame. For the 
r+ — !■ tt+tt'^F mode we exploit the presence of the 
tt" in the final state and the dominance of the decay 
through the p resonance by means of the combined 
quantity Sp: 



^PDG 



,PDG 



rPDG 
P 



-r + i- 



(3) 
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where m^Trf is the reconstructed invariant mass of 
the tt+tt" candidate, the the reconstructed in- 

variant mass of the tt" candidate, m^^^ and T^^^ 
arc the nominal values [6] for the p mass and width, 
m^f^ is the nominal 7r° mass and a^o = 8 MeV is the 
experimental resolution on the 7r° mass determined 
from data. 

We optimize the selection on the B^B~ MC and sig- 
nal MC for the best s/\/s + b, where s is the expected 
signal and b is the expected background from B~^B~ 
events, in the hypothesis of a branching fraction of 
1 X 10~*, for each mode separately. The optimiza- 
tion procedure is performed simultaneously on all the 
discriminating variables in order to take into account 
any correlations. For each discriminating variable ai 
we choose a discrete number A^^ . of possible selection 
criteria. We build N= Na^ x x ....Na^ possible 
selections by combining all the possible cut values for 
different variables. We include E'extra in the optimiza- 
tion, which defines the optimal signal window. We 
choose the selection corresponding to the best value 
of s/V s + b. The optimized selection criteria are re- 
ported in Table II. 



Table II Selection criteria optimized for each t decay 
mode. 



Variable 










i^extra ( GeV) 


< 0.160 


< 0.100 


< 0.230 


< 0.290 


7r° multiplicity 








< 2 


n.a. 


Track multipl. 


1 


1 


< 2 


1 


\cOsdTB 1 


< 0.9 


< 0.9 


< 0.7 


< 0.7 




< 1.25 


< 1.85 


> 1.5 


n.a. 




< 0.9 


n.a 


< 0.5 


< 0.55 




n.a. 


n.a. 


n.a. 


> 1.5 


p quality 


n.a. 


n.a. 


n.a. 


< 2.0 


E^o (GeV) 


n.a. 


n.a. 


n.a. 


> 0.250 



We compute the efficiency as the ratio of the num- 
ber of signal MC events passing the selection criteria 
and the number of events that have a m-^s peaking tag 
B candidate, in the signal region rriES > 5.27 GeV. 
We evaluate the efficiencies on a signal MC sample 
which is distinct from the sample used in the opti- 
mization procedure. A small cross-feed is present in 
some modes and it is taken into account in the com- 
putation of the total efficiency. 

The total efficiency for each selection is: 



Si 



(4) 



where el is the efficiency of the; selection i for the r 
decay mode j, ndec = 7 is the number of of r decay 
modes and fj are the fractions of the r decay mode 
as estimated from the signal MC sample with a recon- 
structed B. Table III shows the estimated efficiencies. 



Table III Selection efficiency (%) for the r decay modes. 
The last row shows the total efficiency including small 
cross-feed from other r decays (not shown explicitly) and 
weighted by the decay abundance at the tag selection level. 
The errors are statistical only. Note that mode dependent 
-Bextra selection is applied as reported in Table II. 



— > e^i/u 
r"*" p.'^w 
r"*" — > TT^V 

T+ ^ 7r+77°!7 


19.3 ± 1.1 
10.8 ± 0.9 
19.7 ± 1.3 
7.0 ± 0.5 


Total: 


9.8 ± 0.3 



In order to determine the expected number of back- 
ground events in the data, we use the final selected 
data samples with i^extra between and 2.4 GeV. 

We first perform an extended unbinned maximum 
likelihood fit to the tues distribution in the final sam- 
ple, in the i?extra sideband region 0.4 GeV < -Eextra < 
2.4 GeV. We use as probability density function(PDF) 
for the peaking component a Gaussian function joined 
to an exponential tail (Crystal Ball), and as a PDF for 
the; non peaking component a phase space motivated 
threshold function [27] (Argus). From this fit, we de- 
termine a peaking yield Tv^^'sii.data signal shape 
parameters, to be used in later fits. Wo apply the same 
procedure on B^B~ MC events which pass the final 
selection and determine the peaking yield N^]f^'^'~' . 
By fitting toes in the i?cxtra signal region (with the 
Crystal Ball parameters fixed to the values determined 
in the i^extra sideband fits) in the MC sample, we de- 
termine the MC peaking yield N^°^'^'^'^\ We perform 
a similar fit to data in order to extract the amount 
of combinatorial background ricomb, as the integral of 
the Argus shaped component in the toes > 5.27 GeV 
domain. 

We estimate the total background prediction in the 
signal region as: 



»t1ow,MC 



high, data 



pk 



N: 



high.MC + "comb 



(5) 



pk 



Due to the low statistics of the final selected sample an 
alternative procedure for estimating the background, 
which avoids fitting the toes distribution in the fi- 
nal samples, is also performed. We fit toes in the 
preselection samples iSextra sideband to determine the 
combinatorial background yield -/V^q^^'^^*'* ^'^^^^^ and 
in the low i^extra region to determine the amount of 
combinatorial background in the preselection signal 
region n^l^. The preselection samples are defined 
by a looser selection only requiring particle identifica- 
tion for the charged track and a lower bound on the 
CM momentum of 0.6 GeV/c for the pion candidate 
and of 1.0 GcV/c for the 7r+7r° candidate. 

We then fit the toes distribution in the final samples 
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in the i?cxtra sideband to determine the combinatorial 
background yield A^^omb'^''*''- ^e estimate the com- 
binatorial background in the low ii'cxtra 
selected sample Ucomb scaling n^J^^^ as 



in the final 



prcscl 
n , X 
comb 



high, data 



high, data prcscl 



N. 



(6) 



The same procedure is applied to MC to subtract 
the combinatorial background and determine the MC 
peaking yield as: 



N. 



low,MC 



pk 



low,MC 



MC prcscl 
'"'comb ^ 



yyhighjMC 
comb 



N, 



high,MC prcscl 



coinh 



(7) 

This peaking component is scaled by the same scale 
factor 7Vphigh,data/^Wgh,MC ^^^^ meihoA and it 

is added to the combinatorial of Eq. (6) to obtain the 
background predictions as in Eq. (5). 

We use the mean value between the two methods to 
obtain the final estimate. We determine the error on 
the estimate by adding in quadrature the statistical 
error from the first method with the half of the differ- 
ence between the two values to obtain the total error 
on the background prediction. The expected numbers 
of background events in the four reconstructed modes 
are shown in Table IV. 



Table IV The mean value of expected number of back- 
ground events in the signal region, using method A and 
method B, respectively. The last column shows the aver- 
age background prediction used in the likelihood scan to 
extract the signal branching fraction. 





method A method B average 


uV 

+ + 0— 


1.1 ±0.9 1.9 ±1.0 1.5 ± 1.4 
1.9 ±1.0 1.6 ±0.8 1.8 ± 1.0 
6.5 ±2.0 7.1 ±1.9 6.8 ± 2.1 
4.1 ±1.4 4.3 ±1.3 4.2 ± 1.4 


All modes 


13.6 ±2.8 14.9 ±3.7 14.3 ± 3.0 



3.2. Results 

After finalizing the signal selection criteria, we mea- 
sure the yield of events in each decay mode in the 
signal region of the on-resonance data. Table V lists 
the number of observed events in on-resonance data 
in the signal region, together with the expected num- 
ber of background events in the signal region. Figs. 7 
and 8 show the E'oxtra distribution for data and ex- 
pected background at the end of the selection. The 
signal MC, scaled to a branching fraction B{B^ — > 
t'^Vt) = 10^"^ is overlaid for comparison. The i^cxtra 
distribution is also plotted separately for each r decay 
mode. 

We determine the 5+ — s- t^Vt branching fraction 
from the number of fitted signal candidates s in data 



Table V The observed number of on-resonance data events 
in the signal region is shown, together with the mean num- 
ber of expected background events. 



r decay mode 


background Observed 


— > e^vV 
r+ ^ ^+^"77 


1.47 ± 1.37 4 

1.78 ± 0.97 5 

6.79 ± 2.11 10 
4.23 ± 1.39 5 


All modes 


14.27 ± 3.03 24 



according to: 



(8) 



where N'^^^ is the number of tag B^ meson correctly 
reconstructed, e^^, £*fg are the tag B efficiencies in 
generic BB and signal events respectively, and Sn, is 
the efficiency to select a signal B ^ tv decay in a 

tagged event. The ratio r^ = was determined 

from MC simulation to be r^ = 0.939 ± 0.007. 

The results from each of our four signal decay modes 



Data 

^^^^ total background prediotion 
^^^^ combinatorial background 

Signal MC (scaled to BR-i ^) 




Figure 7: The iSextra distribution for all the r modes com- 
bined. 





Figure 8: The Sextra 
(b) r+ ^ fi^vV, (c) r 
modes. 



distribution for the (a) r"*" — * e'^vv. 



■R^u, and (d) r+ —* ■k^-k'^v 
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(rich) are combined using the estimator Q = C{s + 
b)/£{b), where £{s + b) and £{b) are the hkehhood 
functions for signal plus background and background- 
only hypotheses, respectively: 

i=l i=l 

(9) 

where Si = £iS. We include the systematic uncer- 
tainties, including those of a statistical nature, on the 
expected background {bi) in the likelihood definition 
by convolving it with a Gaussian function. The mean 
of the Gaussian is b^, and the standard deviation {(Jhi) 
of the Gaussian is the error on bi [28]. 

We compute the central value of the branching 
fraction (including statistical uncertainty and uncer- 
tainty from the background) by scanning over sig- 
nal branching fraction hypotheses between 0.0 and 
3.0 X 10"'' in steps of 0.025 x 10^'' and computing 
the value of C{s + b)/C{b) for each hypothesis. The 
branching fraction is the hypothesis which minimizes 
Q = — 21n(£(s + b)/C{b)), and the statistical uncer- 
tainty is determined by finding the points on the likeli- 
hood scan that occur at one unit above the minimum. 
The systematic error is computed for the branching 
fraction as a fraction of the central value. 

The main sources of uncertainty in the determina- 
tion of the B'^ —> t'^Vt branching fraction come from 
the estimation of the tag yield and efficiency, the re- 
construction efficiency of the signal modes and the 
number of expected background events. We estimate 
a systematic uncertainty due to the teclmiriue used 
to determine the tag B yield and reconstruction ef- 
ficiency, varying the MC based criteria to define the 
B^B^ non-peaking component of the jties shape. Ob- 
serving the effect due to the different choices on the 
branching fraction measurement, we assign a system- 
atic uncertainty of 3%. The systematic uncertainty 
of the signal efficiencies depends on the r decay mode 
and includes the effects of the tracking of charged par- 
ticles, particle identification, and the 7r° reconstruc- 
tion efficiency. The dominant efficiency uncertainty 
was found to be due to Eextra selection requirement 
and is found to be 15%. The contributions due to the 
tracking and to the choice of the Eextra selection are 
treated as corrclatcid among the different channels. 

The systematic uncertainty on the background pre- 
diction has been estimated using two different meth- 
ods to extract the combinatorial background, and it 
has been found to be negligible with respect to the 
statistical uncertainty. This uncertainty has been in- 
corporated in the likelihood definition used to extract 
the branching fraction. 

We determine the branching fraction central value 
to be 

B{B+ T+Ur) = 1.8lo^(stat-hbkg)±0.3(syst))xlO-'*. 

(10) 



We obtain a significance of 2.7 a from S = y^Q,„i„, 
where Qmjn is the minimum value of the likelihood 
ratio, if we do not include the uncertainties on the 
background predictions. If we include the background 
uncertainty into account we obtain the smaller signifi- 
cance of 2.2 a. We compute the 90% C.L. upper limit 
using the CL^ method [29] to be B{B+ r+Vr) < 
3.4 X 10-4. 

3.3. Combined Results 

We report here the branching fraction obtained by 
combining the hadronic tag analysis result, described 
in this document, with the BaBAR semileptonic tag 
analysis result, based on a statistically independent 
data sample and reported in [30]. Both the analy- 
ses proceed in a similar manner for the four r decay 
mode reconstruction and use the same likelihood ra- 
tio scan technique to combine the different t decay 
modes. Fig. 9 shows the ii'extra distribution for all the 
T modes combined in the analysis with semileptonic 
tag. 




ri F T"!''''?;',':',^','!^^ 

H) 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 9: The Ee^tia. distribution for all the r modes 
combined in the analysis with semileptonic tag. The 
dashed line shows predicted signal distribution assuming 
the — » r^!/ branching fraction of 10"*. 



The likelihood ratio is extended to combine eight 
reconstruction modes, four from the hadronic tag and 
four from the semileptonic tag. We obtain 

B{B+ T+Vr) = 

(1.20tO;*°(stat)toj^(bkg)±o.22(syst)) x 10-^ (11) 

where the first error is statistical, the second is due 
to the background prediction systematic uncertainty, 
and the third one is due to the other systematic 
sources. We obtain a significance of 2.6 a including 
the uncertainty on the expected background (3.2 cr if 
this uncertainty is not included). 
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Fig. 10 shows the C{s + b)/£{b) scan as a function of 
signal branching fraction (the background uncertainty 
is included in the likelihood). 



^J5 




Hadronic 

Semileptonic 
Combined 



xlO" 



0.05 O.L 0.15 0.2 0.25 0.3 0.35 0.4 

Branching Fraction 



Figure 10: Scan of Q — ~2ln{jC{s + b) / £{b)) as a function 
of the — > T^^i/ decay branching fraction. The uncer- 
tainty on the expected number of background events has 
been included in the likelihood definition. The results are 
shown for the hadronic tag, semileptonic tag, and com- 
bined analysis. 



3.4. Discussion 

In summary, we report here preliminary branching 
fraction of the t^v decay using two inde- 

pendent methods of tag B reconstruction: hadronic 
and semileptonic B decay in the recoil. We find 
a 2.6 tr (3.2 a not including expected background 
uncertainty) excess in data which can be converted 
to a preliminary branching fraction central value of 
B(B± ^ ^±^) = (1.20t°:«t°:fo ± 0.22) X 10-4. 

Given the measurement of the B^ t^v branch- 
ing fraction, we may constraint CKM parameter Vub 
using Eq. (1). See Fig 11 for illustration. Alterna- 
tively, given the best world measurements of param- 
eters present in Eq. (1), we can constraint tan/3 and 
charged Higgs mass in the model with two Higgs dou- 
blets, as shown in Fig 12. 



0.5 



-0.5 



-1.5 




Constraint from B* — ^ t* and Am^ 

I I I I I I I I I I rrr , , , , ittt 



-0.5 



p 



1.5 



Figure 11: Constraints at 95% C.L. on p and r] parameters 
(through measurement of Vuh) of the CKM triangle using 
T^u results presented here, see Eq. (1). Constraints 
from Bd-mixing are also shown, along with combined re- 
sults [31]. 

.T- 500 



> 450 
S 400 
3^350 
300 
250- 
200 
150 
100 r 




10 20 30 40 50 60 70 80 90 100 

tanp 

Figure 12: Constraints on tan/3 and charged Higgs mass 
from B^ T^u results presented here, see Eqs. (1) and 
(2). Shaded areas are excluded at 95% C.L. 



4.1. Analysis Method 



4. B± v?(1020)K*(892)± 

One of our goals is to test the quark-spin projections 
in the B (pK* decay, as shown in Fig. 3. However, 
equivalently we can measure spin projections of the 
two vector mesons in the B decay. There are three 
complex amplitudes for the three spin projections pos- 
sible, and therefore there are a total of 12 non-trivial 
independent real parameters describing the B and B 
decays. Information about the amplitude sizes and 
phases can be deduced from analysis of angular dis- 
tribution as illustrated in Fig. 13. 



The B± -> Lp{1020)K*^ -> {K+K^){Kn)^ can- 
didates are analyzed with two (Kn)^ final states, 
KgTT^ and K^tt^. The neutral pseudoscalar mesons 
are reconstructed in the final states Kg tt+tt" and 
^ ry^^ define the hclicity angle 9i as the an- 
gle between the direction of the K or K'^ meson from 
K* Kit iOi) or ip ^ K+K- (62) and the direc- 
tion opposite the B in the K* or (p rest frame, and 
$ as the angle between the decay planes of the two 
systems, see Fig. 13. The differential decay width has 
four complex amplitudes Ajx which describe two spin 
states of the Ktt system (J = 1 or 0) and the three 
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Figure 13: Definition of decay angles given in the rest 
frames of the decaying parents. 



helicity states of the J = 1 state (A = or ±1): 



dUidn2d<^ 



^AjAl'j(Wi,$)rf\-7i2,0) 



^2) 



where Hi = cos6'i and Yj are the spherical har- 
monics with J = 1 for ic*(892) and J = for 
{K'it)q. We reparameterize the amplitudes as Ai±i ~ 
(Ai||±^i^)/\/2. 

We identify B meson candidates us- 
ing two kinematic variables: "^es = 
[(s/2 + pr-PB)VS|-p|]'/' and AS 
[EtEb - Pr • Pb - s/2)/Vs, where {Eb,Vb) 
is the four-momentum of the B candidate, and 
(i?r,Pr) is the e~'"e~ initial state four-momentum, 
both in the laboratory frame. We require tties > 5.25 
GeV and jASj < 0.1 GeV. The requirements on 
the invariant masses are 0.75 < ruK-K < 1-05 GeV, 
0.99 < TO^^ < 1.05 GeV, |to^^ - m^o| < 12 MeV, 



< 150 MeV for the K*^, ip, K^, Vi- 



and 120 < 

and 7r°, respectively. For the candidates, we also 
require the cosine of the angle between the flight 
direction from the interaction point and momentum 
direction to be greater than 0.995 and the mea- 
sured proper decay time greater than five times its 
uncertainty. 

To reject the dominant e+e^ — > quark-antiquark 
background, we use the angle 9t between the B- 
candidate thrust axis and that of the rest of the event, 
and a Fisher discriminant T [32]. Both variables are 
calculated in the center-of-mass frame. The discrim- 
inant combines the polar angles of the S-momentum 
vector and the B-candidate thrust axis with respect 
to the beam axis, and two moments of the energy flow 
around the B-candidate thrust axis [32]. 

To reduce combinatorial background with low- 
momentum tt'^ candidates, we require Tii < 0.6. 
When more than one candidate is reconstructed, 
which happens in 7% of events with Kg and 17% with 
tt", we select the one whose of the charged-track 
vertex fit combined with of the invariant mass con- 
sistency of the Kg or tt^ candidate, is the lowest. We 



define the 5-quark flavor sign Q to be opposite to the 
charge of the B meson candidate. 

We use an unbinned, extended maximum-likelihood 
fit [11, 15] to extract the event yields Uj and the pa- 
rameters of the probability density function (PDF) 
Vj . The index j represents three event categories used 
in our data model: the signal B^ (p{Kn)^ {j — 1), 
a possible background from B^ — > /o(980)ii'*^ {j = 
2), and combinatorial background (j = 3). The super- 
script k corresponds to the value of Q = ± and allows 
for a CP-violating difference between the B'^ and B~ 
decay amplitudes {A and A). In the signal category, 
the yield and asymmetry of the B^ LpK*{S'd2)^ 
mode, Usig and Acp, and those of the B^ (^{K'k)^^ 
mode are parameterized by applying the fraction of 
LpK*{892)^ yield, iJ' , to n\. Hence, n^ig = n'^ x + 
X /i^, Acp = [n^ x — x ii~)/nsig, and the 
iy9(i^7r)Q^ yield is rii x (1 — /i+) + x (l — ^^). 

The likelihood Ci for each candidate i is defined 
as A = Ej,^^^ M^ C, I), where the PDF is 

formed based on the following set of observables = 
{Til, H2, ruK-K, rrij^j^, AE, toes, ^, Q} and the 
dependence on and polarization parameters C is 
relevant only for the signal PDF Vi- The remaining 
PDF parameters ^ are left free to vary in the fit for the 
combinatorial background and are fixed to the values 
extracted from Monte Carlo (MC) simulation [23] and 
calibration B — > Dn decays for event categories j — 1 
and 2. 

The helicity part of the signal PDF is the ideal an- 
gular distribution from Eq. (12), multiplied by an em- 
pirical acceptance function 5(7Yi,7Y2,$) = GiCHi) x 
02(^2) ■ Here, the amplitudes Aj\ are expressed in 
terms of the polarization parameters C = {II, </'||) 
^0, A%, Acp, A^||, A(j)±, A60} defined in Ta- 



ble VI. CP-violating differences are incorporated via 
the replacements in Eq. (12) for B~^ decays: /l 
fL X (1 +^o,p X Q), f±^f±x {1+A^p X Q), ^ 
((/.|| + A(/)|| X Q), -> (0i + tt/2 + (A0_L + 7r/2) x Q) , 
and i5o {So + A60 x Q). 

A relativistic spin-J Breit-Wigner amplitude pa- 
rameterization is used for the resonance masses [6, 33] , 
and the {KTr)^^ rnx-K amplitude is parameterized 
with the LASS function [20]. The latter includes 
the /-Cq (1430)^ resonance together with a nonresonant 
component. The interference between the J = and 1 
(i^Tr)* contributions is modeled with the three terms 
27?,e(^iA^oo) ^'i- (12) with the four-dimensional 
angular and to^tt parameterization and with depen- 
dence on /i*^ and C- 

The signal PDF for a given candidate i is a joint 
PDF for the helicity angles and resonance mass as 
discussed above, and the product of the PDFs for 
each of the remaining variables. The combinatorial 
background PDF is the product of the PDFs for inde- 
pendent variables and is found to describe well both 
the dominant quark-antiquark background and the 
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Table VI Summary of results for the ipK* (892)^ decay. The twelve primary results arc presented for the two 

decay subchannels along with the combined results, where the branching fraction B is computed using the number 
of signal events risig and the total selection efficiency e, which includes the daughter branching fractions [6] and the 
reconstruction efficiency Ereco obtained from MC simulation. The definition of the six CP-violating parameters allows 
for differences between the S"*" and B~ decay amplitudes A and A with superscript Q = — and +, respectively. The 
systematic uncertainties are quoted last and are not included for the intermediate primary results in each subchannel. 
The dominant fit correlation coefficients (C) are presented, where we show correlations of So with <p\\/(j)± and of ASo with 
A0||/A0x. 



parameter 


definition 


B 


r/Ftotal 


h 


|Aio|VS|^iA 


f± 


|>1ix|VE|^ia 


h 


arg(^i|l/^i„) - 


4>± — n 


arg(AixMio) - 


5o —TT 


arg(^ooMio) - 



K*{892y 



A'*(892)± ^ 



combined 



- 7r 

TT 



(10.5 ± 1.4) X 10" 
0.51 ±0.07 
n 22+°-°'^ 

"•^^-0.06 

+0.28 



-0-75 
-0.15 ±0.24 
-0.25 ± 0.24 



(11.6 ± 1.5) X 10" 



-0.09 
I +0.09 
-0.08 



o.2i: 

-0.77 ± 0.35 

-0.89l°;f6 
-1-0.11 ±0.31 



(11.2 ±1.0 ±0.9) X 10"*' 

0.49 ± 0.05 ± 0.03 

0.21 ± 0.05 ± 0.02 
-0.67 ±0.20 ±0.07 
-0.45 ± 0.20 ± 0.03 
-0.07 ±0.18 ±0.06 ±37%/ ±36% 



}-58% 
} ±56% 



Acp (r+ -r-)/(r+±r-) 
ift-fD/Ut + fl) 

Ahp {ft - fl)/(.ft + fl) 

A0II (<^+-0-)/2 
A0X (^+ 7r)/2 



Z 

^reco 



-0.09 ±0.13 
±0.24 ± 0.15 
±0.12 ±0.31 
±0.02 ±0.28 
±0.18 ±0.24 
±0.13 ± 0.24 

102 ± 13 ± 6 
(2.53 ±0.13) % 
(22.3 ± 1.2) % 



±0.07 ±0.13 
±0.09 ± 0.20 

in 41+0-54 

±0.22 ±0.35 

±0.48lS:fo 
±0.34 ± 0.31 

11711^ ± 7 
(2.59 ±0.17) % 
(16.0 ± 1.0) % 



0.00 ± 0.09 ± 0.04 
±0.17 ±0.11 ±0.02 
±0.22 ± 0.24 ± 0.08 
±0.07 ± 0.20 ± 0.05 
±0.19 ±0.20 ±0.07 
±0.20 ±0.18 ±0.03 



}-50% 

} ±57% 
±37%/ ± 37% 




-0.1 0.1 5.25 5.26 5.27 5.28 5.29 



AE (GeV) nigs (GeV) 

Figure 14: Projections onto the variables (a) rxiK-n, (b) 
rrij^j^, (c) AS, and (d) rriES for the signal B^ — > (p^Kn)"^ 
candidates with a requirement discussed in the text. 
The solid (dashed) lines show the signal-plus-background 
(background) PDF projections. 



background from random combinations of B tracks. 
The signal and background PDFs are illustrated in 
Figs. 14 and 15. For illustration, the signal frac- 
tion is enhanced with a requirement on the signal- 
to-background probability ratio, calculated with the 




Q <I> (degrees) Q <t> (degrees) 

Figure 15: Projections onto the variables (a) Ti.i, (b) Ti.2, 
(c) Q and (d) the differences between the Q $ projec- 
tions for events with Tli H2 > and with Hi TI2 < 
for the signal B^ — » ip{KiT)^ candidates following the 
solid (dashed) line definitions in Fig. 14. The step in the 
"Hi PDF distributions is due to the selection requirement 
Hi < 0.6 in the ipiK'^Tr") channel. 



plotted variable excluded, that is at least 50% effi- 
cient for signal B^ ip{KTr)^ events. We use a sum 
of Gaussian functions for the parameterization of the 
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2 4 6 

^11 (rad) 



Figure 16; Distribution of 2 In £ as a function of for the 
two fit configurations: including P- and S-wave Kn inter- 
ference into account (solid points) and ignoring it (open 
points). The former resolves the phase ambiguity. 




2 4 6 

O^(rad) 



Figure 17: Distribution of 2 In £ as a function of (j)± for 
the two fit configurations: including P- and S-wave Ktt in- 
terference into account (solid points) and ignoring it (open 
points). The former resolves the phase ambiguity. 



ated systematic uncertainties. We allow for a flavor- 
dependent acceptance function and reconstruction ef- 
ficiency in the study of asymmetries. The biases from 
the finite resolution of the angle measurements, the 
dilution due to the presence of fake combinations, or 
other imperfections in the signal PDF model are esti- 
mated with MC simulation. 

The nonresonant K'^ K~ contribution under the 
is accounted for with the JqK*^ category. 

Its yield is consistent with zero. The ruj^j^ PDF 
shape in this category is varied from the resonant 
to phase-space and the yield is varied from the ob- 
served value to the extrapolation from the neutral 
i?-decay mode [15] to estimate the systematic uncer- 
tainties. Additional systematic uncertainty originates 
from other potential B backgrounds, which we esti- 
mate can contribute at most a few events to the signal 
component. The systematic uncertainties in efficien- 
cies are dominated by those in particle identification, 
track finding, and and 7r° selection. Other sys- 
tematic effects arise from event-selection criteria, 
and K*^ branching fractions, and the number of B 
mesons. 

4.3. Discussion 

The yield of the ip{KiT)'^ contribution is 57l}3 
events with a statistical significance of 7.9cr, com- 
bining the |AooP term and the interference terms 
2R.e{AixAQ^) , which confirms the significant S'-wave 
Kt: contribution observed in the neutral -B-decay 
mode [15]. The dependence of the interference on the 
Ktt invariant mass [15, 20, 21] allows us to reject the 
other solution near (27r — ^||,7r — relative to that in 



signal PDFs for /S.E, mES, and J-. For the combinato- 
rial background, we use polynomials, except for tties 
and J- distributions which are parameterized by an 
empirical phase-space function and by Gaussian func- 
tions, respectively. Resonance production occurs in 
the background and is taken into account in the PDF. 

4.2. Results 

We observe a nonzero B^ ipK*{d,^2)^ yield 
with significance, including systematic uncertainties, 
of more than 10a. The significance is defined as the 
square root of the change in 2 In £ when the yield is 
constrained to zero in the likelihood C. In Table VI, 
results of the fit are presented, where the combined 
results are obtained from the simultaneous fit to the 
two decay subchannels. 

We repeat the fit by varying the fixed parameters 
in ^ within their uncertainties and obtain the associ- 




0.75 1.45 2.15 

mK,(GeV) 



Figure 18: Distribution of the Kn invariant mass extended 
above 1.6 GeV from the study of B° (f){K+Tx-) de- 
cays in Ref. [15]. The data distribution is shown with 
a requirement to enhance the signal as discussed in re- 
gard to Fig. 14 in text. The solid (dashed) line shows 
the signal-plus-background (background) expected distri- 
butions. The arrows indicate the higher mass range, 1.60 
to 2.15 GeV. 
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Table VII Summary of branching fraction (B) and longitudinal polarization fraction (/l) results for the various B —* tpK* 
decay on BABAR,. The quantum numbers are given for the K* resonances. For a complete list of up to 12 independent 
parameters measured in the (pA"*(892)'' and iy5i('|(1430)'' decay modes see Ref. [15]. 





B decay mode 


Branching (10"'') 


fL 


0+ 

1" 
1" 
1" 

2+ 
3" 

4+ 


(^^-0(1430)" 
(/pit"* (892)" 
^K*{m2)+ 

<^ ATI (2045)" 


4.6 ±0.7 ±0.6 
9.2 ±0.7 ±0.6 
11.2 ±1.0 ±0.9 
<3.5 (0.71J:? ± 1.1) 
7.8 ±1.1 ±0.6 
< 2.7 (-0.9 ±1.4 ±1.1) 
< 15.3 (6.0j:to ±4-1) 


0.51 ±0.04 ±0.02 
0.49 ±0.05 ±0.03 

0.85to;;;? ± 0.04 




Figure 19: Feynman diagram describing the B tpK* de- 
cay with the new mechanisms within the standard model, 
such as annihilation penguin [18] (top) or QCD rescatter- 
ing [19] (bottom). 



Table VI with significance of 6.3cr, including system- 
atic uncertainties. Figs. 16 and 17 illustrate resolution 
of the phase ambiguity. 

The iy — A) structure of the weak interactions, he- 
licity conservation in strong interactions, and the s- 
quark spin flip suppression in the penguin decay di- 
agram suggest \Aiq\ > |^i+i| > [10]. This 
expectation disagrees with our observed value of /l. 
We obtain the solution </)|| ~ without discrete am- 
biguities, which is consistent with the approximate 
decay amplitude hierarchy \Aiq \ ~ |^i+i| ^ 

We find that and 0|| deviate from either tt or 
zero by more than S.lcr and 2.4(t, respectively, in- 
cluding systematic uncertainties. This indicates the 



presence of final-state interactions not accounted for 
in naive factorization. Our measurements of the six 
CP-violating parameters are consistent with zero and 
exclude a significant part of the physical region. We 
find no evidence of CP violation in this decay. 

We have performed a full amplitude analysis and 
searched for CP-violation in the angular distribution 
of the (pK*^ decay. Our results are summa- 

rized in Table VI. Using similar techniques, we have 
also searched for the (j)K*^ decays with the 

tensor Kl[n%Qf and Kl{2Mbf , vector if* (1680)°, 
and scalar nonresonant {Ki:)'^' contributions with 
A'*" A'+TT" invariant mass above 1.6 GeV. Our 
results are summarized in Table VII. Fig. 18 shows 
distribution of the ATtt invariant mass in this study. 
We do not find significant signal with the above reso- 
nances and place upper limits on the _B-decay branch- 
ing fractions. 

These results find substantial Ai^i amplitude in the 
B ipK*{892) decay and point to physics outside the 
standard model as shown in Fig. 4 [17], where scalar 
interaction would look more consistent with the data, 
or new dynamics as shown in Fig. 19 [18, 19], where 
penguin annihilation is more consistent with sizable 
Ai+i and small Any satisfactory solution to 

this polarization puzzle is expected to explain all po- 
larization data and to be predictive for other experi- 
mental measurements. 
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